Instabilities of waves in nonlinear disordered media 
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We develop a self-consistent theory of temporal fluctuations of a speckle 
pattern resulting from the multiple scattering of a coherent wave in a weakly 
nonlinear disordered medium. The speckle pattern is shown to become un- 
stable if the nonlinearity exceeds a threshold value. The instability is due to 
a feedback provided by the multiple scattering and manifests itself in sponta- 
neous fluctuations of the scattered intensity. The development of instability 
is independent of the sign of nonlinearity. 
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Interplay between nonlinearity and multiple scattering has received considerable atten- 
tion during, the recent years. Second harmonic generation pi phase conjugation,! and optical 
Kerr effecta have been studied in the weak scattering limit {k^ ^ 1, — the wavenumber 
in the medium, £ — the mean free path). On the other hand, the nonlinearity is shown to 
modifv the coherent backscattering peak,l and to affect the angular correlations of scattered 
waves.a 

In the present Letter, we consider the temporal fluctuations of waves in weakly nonlin- 
ear disordered media. We show that, in contrast to linear media, where the fluctuations of 
scattered waves with time are due to scatterer motionJiHl the wave scattered in a nonlinear 
medium may exhibit temporal fluctuations even if the scatterers are immobile. These fluc- 
tuations are spontaneous and represent a signature of instability of the multiple-scattering 
speckle pattern. Although unstable regimes are rather common for nonlinear waves in feed- 
back systems,EHli our study is the first to consider the (multiple) scattering as a feedback- 
providing mechanism in a nonlinear medium. The long-range spatial correlations of intensity 
fluctuations appear to be crucial for the feedback efficiency. 

We consider a scalar plane wave of frequency uj incident upon a nonlinear disordered 
medium occupying the half-space 2 > 0. Wave amplitude in the medium V'(r, t) obeys a 
nonlinear wave equation: 

{ + kl [eo + Seiv) + |V(r, t) \^] } V(r, t) = J(r, t), (1) 

where /cq is the free-space wavenumber, Eq = Uq is the average dielectric function, Se{r) is the 
fluctuating part of the dielectric function, 82 is a nonlinear constant, and J{r,t) is a (slowly 
varying) source term [\{d/dt)J{r,t)\ <C a;|J(r, t)|]. Eq. (|I]) describes propagation of waves 
in centrosymmetric media, where nonlinear terms quadratic in are forbidden.0 Equations 
of the same type appear for electromagnetic waves in plasma,Ej as well as in the context of 
the theory of hydrodynamic turbulence .0 A nonlinear Schrodinger equation similar to Eq. 
(0) has been recently considered by Spivak and Zyuzin.El 
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We assume the source term J(r, t) in Eq. ([^) to consist of two parts: Jo{r), representing 
an ideal plane wave incident at z = 0, and 6J{r,t), which is a random noise term. In the 
present paper, we are only interested in the case of infinitely weak noise: \6J\ — 0. In a linear 
medium (£2 = 0), the solution of Eq. is then unique and time independent: ip{r, t) = ip^r). 
Consequently, the autocorrelation function (7i(r,r) = {ip{r,t)i/j*{r,t + t)) / (/(r, t)) = 1, 
where /(r,t) = \tMr,t)f is the wave intensity. If 62 7^ 0, however, Eq. (0) may admit 



multiple solutions,li3 and (even being infinitely weak) the noise 6 J can provoke the system 
to jump from one solution to another. Obviously, such a behavior will result in spontaneous 
fluctuations of ip{r,t) with time, yielding gi{r,T) < 1. In the general case, the calculation 
of gi appears to be too involved. We therefore limit ourselves to long times r ^ Tj^mp, 
where Tjump is the typical time needed for the system to jump from one solution to another. 
Calculation of Tjump is beyond the scope of our analysis, but we expect Tjump to be inversely 
proportional to the speed of wave in the medium. For r 3> Tj^mp, should become 

independent of r, and its asymptotic value can be considered as a measure of the mutual 
resemblance of solutions of Eq. ([^). To simplify the further analysis, we adopt the following 
four assumptions: 

(i). Propagation of waves in a weakly nonlinear disordered medium is diffusive with a mean 
free path i unaffected by nonlinearity. This implies that nonlinear refraction is negligible at 
distances of order £, and consequently, that An^ki <^ 1, where An = n2/o, n2 = e2/{2nQ), 
and Jo = (/(r, t)) in the absence of absorption. 

(m). Statistics of the wave field scattered in a weakly nonlinear disordered medium 
is close to Gaussian. Consequently, the factorization approximation (72(1", t) = 
(/(r,t)/(r,t + r))/(/(r,t))^ ~ 1 + | 5fi(r, r)|^ applies.! This assumption is justified under 
the same conditions as (z), since the Gaussian statistics of the total scattered wave field 
ip{r, t) is a consequence of the complete randomization of phases of partial waves arriving at 
r. The reason for the randomization is that the typical distance i between individual scatter- 
ing events in a multiple-scattering sequence is much larger than the wavelength {k£ ^ 
Obviously, such a mechanism of phase randomization is equally effective in both linear and 
weakly nonlinear media, as long as (i) holds. 

{in). Intensity of the third harmonic remains always much smaller than the intensity of the 
fundamental wave. This implies that the medium has a sufficient degree of dispersion for 
the phase matching condition^ to be violated: |AA;| £ ^ 1 with Ak = k^ — 3k and k^ being 
the wavenumber at frequency 3uj. 

(iv). Functional form of the r-dependence of gi is the same as in the linear medium: 
^i(r,r) = exp[— /5(r) z/£], where /5 <C 1 is some unknown function. To justify this ansatz, 
we make use of (i) and write gi as a sum of contributions of paths of various lengths s:ili 



Here P{z,s) oc zs~^^'^ ex'p[—3z'^ / (Ms)] is the relative "weight" of paths of length s, and 
(A</)^(r))^ is the variance of the phase difference between two waves propagating along the 
same path but separated in time by r, averaged over all possible realizations of disorder 
and over all paths of the same length s. From here on, we denote such an averaging by 





dsP(z,s). 



(2) 
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(■ ■ ■)^. As a function of s, P{z, s) is peaked at sq = /{2t). The second equality of Eq. 
(@) is obtained by assuming that (A(^^(r))^ increases monotonically with s and becomes 
of order unity at s = Si(r) ^ sq. The result of integration in Eq. can be written as 
1 - I3{t) zjt ~ exp[-/?(r) zl^\ with ^zji ~ ^(^Si)-^/^ ^ (2so/si)i/2 ^ i_ ^j^g opposite 
limit of ^zji ^ 1, (7i(r, r) vanishes and the functional form of its z-dependence is of no 
importance. Anyway, it will be seen from the following that the exact functional form of the 
r-dependence of (7i(r, r) is not of crucial importance, since g\ is integrated over the whole 
medium during the calculation of the correlation function of diffusely reflected wave. 

Consider now a typical diffusely reflected wave path of length s. The phase </?(t) accu- 
mulated along such a path is 

V?(t) = fco /" ds' {no + n2l[r(s'), t]} , (3) 

where the integration is along the path. The variance of the phase difference A(y9(r) = 
+ t) — ^{t\ averaged over all possible paths of the same length s, is found directly from 
Eq. as a sum of two contributions: 

(Av.^(r)>f ^ = (vr/no) k.tnl {Aliv, r)^)^ s/i, (4) 
(A^2(r)>f = (koirnl (A/(ri, r)A/(r2, r)), (5) 

where AJ(r, r) = /(r, t + r) — /(r, t). {Aip"^ {t))^^^ and {Aip'^{T))f^ originate from the short- 
and the long-range spatial correlations of A/, respectively. 

If r ^ Tjump, the calculation of the r.h.s. averages in Eqs. (||) and is straightfor- 
ward. To evaluate (AJ(r, r)^)^, we integrate (AJ(r, r)^) = 2(J(r, — gf{r,T)] over 
z > with a weight function Ps(r) describing the density of paths of length s in the 
medium.! To find the average of Eq. (^), we calculate the long-range correlation function 
(A/(ri, r)AJ(r2, r)) using the Langevin approachS with the time correlation function of 
Langevin random sources being proportional to gf{r, r). The result is then averaged over all 
pairs of points ri, r2 belonging to the same path with a weight function ps(ri,r2) which is 
calculated similarly to Ps(r). As gi = exp(— enters into the derivations of (A/(r, r)^)^ 
and (A/(ri, r)A/(r2, r))^, {Aip'^^r))^^^'^'^^ appear to depend explicitly on /3. The details of 
the calculations will be presented elsewhere, the final results are: 

(Ay,2(r)>f) = 2n koiAn'H (/^^^ ) j, (6) 



[A^\T)f^ ^ QAn^ X I ^ ['Jf^ ^vW^ ^ 1' (7) 



where H{x) = ^/7^x exp{x'^)[l — Erf(x)], and rio = 1 is assumed for simplicity. Applying Eq. 
we obtain the autocorrelation function of diffusely reflected wave: 



POO 

^i(£,r)= exp[-/j(r)] = / dsP(£,s)exp 

Jo 



F[/3(r)], (8) 



where (A(y9^(r))^ is a sum of the terms given by Eqs. (|^) and (|^). Since the integral of Eq. 
(0) does not contain r in explicit form, we conclude that jS and gi are independent of r. 
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as expected for r ^ Tjump- From here on, gi without arguments denotes the value of the 
autocorrelation function of diffusely reflected wave at r 3> Tj^mp- 

In Eq. (|^), we have introduced a new function F{(3) which depends on An and k^l as 
well. The underlined part of Eq. (§) can be considered as a nonlinear equation for /5, or, 
equivalently, for gi since /3 = —In (71 by definition. Making a simultaneous plot of F{(3) 
and exp(— /3) (see the lower solid curve in the inset of Fig. |lD, we find that Eq. (|]) has two 
solutions: /3 = and (3 > 1. The first solution is "unstable" with respect to the weak noise 
5 J in Eq. (P, since 5 J 7^ introduces an additional dephasing term into Eq. (H), yielding 
F(0) < 1. It is then the second solution (3 > 1 which is the physically realizable one. As 
/? > implies gi < 1, our results show that ipiv.t) exhibits spontaneous fluctuations, and 
that the speckle pattern resulting from the multiple scattering of coherent wave in a weakly 
nonlinear disordered medium is unstable. 

The value of gi, characterizing the magnitude of spontaneous fluctuations {gi = 1 in the 
absence of fluctuations), can be found by a numerical solution of Eq. (see the dashed 
line in Fig. [l|), as well as analytically. Assuming 1 — (71 ^ 1, and noting that the ultimate 
contribution to (A(/)^(r))g is given either by Eq. (|^) [for s/i (^o^)^] or by Eq. (|^ [for 
s/i > (fco^)^], we obtain 
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Up to now, we have considered a nonabsorbing medium. Absorption, however, is in- 
evitable in a real medium and should therefore be necessarily taken into account. This can 
be done by allowing Eq in Eq. (ffl) to have a small imaginary part. Eq. (0) then becomes 



(Av.^(r)> 



(1) 



27r kniAn^ { H 
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where La = [^/(Sfcolmeo)]^''^ ^ £ is the absorption length, and ng = Re^o = 1- Eq. (^ is 
still approximately valid for {i/ La){s/ dy/"^ < 1, while for (£/La)(s/^)^/^ > 1 we get 



,(2) 



6An^ X 



(La/if, > 1. 



(11) 



Graphical solution of Eq. (|[) with P{i,s) oc exp[— 3£/(4s) — si/{3Ll)] is presented 
in the inset of Fig. |^ for several absorption lengths La . As follows from the figure, the ab- 
sorption has to be weak for the second solution of Eq. (||), /5 > 0, to exist. Sufficiently strong 
absorption prevents spontaneous fluctuations of V'(r,^), stabilizing the speckle pattern. The 
condition for the onset of the speckle pattern instability is {d / d(3)F{j3) = —1 evaluated at 
P = 0. This yields: 



p = An^ 



ko£ + 



1, 



(12) 



where we introduce a control parameter p. U p < 1, gi = 1 and the multiple-scattering 
speckle pattern is stable, while for p > 1, an instability shows up leading to (71 < 1 (see 
Fig. IID. While gi characterizes the mutual resemblance of solutions of Eq. (|I|), the number 
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of solutions grows exponentially with p^/^, as can be estimated using the method of Ref. 
|T3|. The striking feature of Eq. (^) is that p can reach unity even for very small |An| 
provided that the extensive parameter La/i is large enough. The stability condition p < 1 
is therefore never satisfied in a nonabsorbing medium. Eq. (|T2]) is consistent with the result 
of Spivak and Zyuzinjll who have shown that the perturbation analysis of the sensitivity of 
speckle pattern in a nonlinear disordered medium to changes of scattering potential fails for 
An'^{L/i)^ > 1, where L is the typical size of the medium. 

The physical origin of the instability of speckle pattern for p > 1 can be revealed by 
realizing that the system "coherent wave + nonlinear disordered medium" has a positive 
distributed feedback. An infinitely weak perturbation of the wave intensity /(r, t) produces a 
change of the local refractive index, altering the phases of waves propagating in the medium, 
and affecting their mutual interference. Since it is this interference which determines I{r,t), 
the feedback loop is built up and closed. For p ^ 1, the feedback is sufficiently strong to 
compensate for the (diffusive on average) spreading of the initial intensity perturbation, and 
the speckle pattern /(r,t) becomes unstable. 

The above analysis can be generalized to the case of Se which varies slowly with time 
(moving scatterers). An additional contribution to (Ay9^(r))^ arises in this case. If, for 
instance, the scatterers undergo Brownian motion with a diffusion coefficient D, this contri- 
bution is 

<A^^(r)>f^ = ^(l + 2An 
^0 L 

where both tq = (4A;qD)~^ and r are assumed to be much larger than Tjump- The time auto- 
correlation functions of diffusely reflected waves obtained by solving Eq. (^ with (A(y9^(r))^ 
given by a sum of contributions of Eqs. (0), ([l^), (O)) and (|13]) are shown in Fig. 0. The 
two lower curves correspond to the unstable regime {p > 1), and therefore (71 (r) < 1 even 
if t/tq = 0. For the upper curve, the speckle pattern is stable {p < 1) and (71 (r) = 1 at 
r/ro = 0. 

In conclusion, we have shown that the multiple-scattering speckle pattern may become 
unstable if the disordered medium exhibits a weak nonlinearity of Kerr type. The insta- 
bility is due to the distributed feedback provided by the scattering, and its development is 
independent of the sign of nonlinearity. The instability threshold decreases as absorption 
is reduced, and the speckle pattern is never stable in a nonabsorbing unbounded medium. 
Finally, we note that our results may be of particular importance in the context of "random 
laser" ,t3 where absorption may be efficiently compensated by stimulated emission. 
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FIG. 1. Main plot: a "bifurcation diagram" of the speckle pattern resulting from the multiple 
scattering of a coherent wave in a semi-infinite nonlinear disordered medium at fixed fco^ = 100. 
The inverse absorption lengths l/La are indicated near each curve, dashed line corresponds to 
a nonabsorbing medium, dotted line shows gi = 1 (linear case). In the absence of absorption, 
the speckle pattern is unstable (171 < 1) at any An 7^ 0, while in a dissipative medium, |An| 
should exceed a certain threshold value for the instability to develop. The threshold values of |An| 
following from Eq. (12) are shown by arrows. Inset: the graphical solution of Eq. (^ for fixed 
|An|, /co^, and the values of l/La indicated near each curve. The dash-dotted line is exp(— /3). 



FIG. 2. The field autocorrelation function of a wave diffusely reflected from a semi-infinite 
nonlinear (solid curves) or linear (dotted curves) disordered medium for fixed An, UqI, and inverse 
absorption lengths I /La indicated near each curve (Tjurnp/To — > is assumed). For the two largest 
values of absorption length, the "nonlinear" speckle pattern is unstable {gi < 1 at t/tq = 0). 
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